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I. I NT 50 DUCT IPS 



An aircraft versus missile engagement is traditionally 
described in two phases, the missile flyout and the Endgame. 
That portion of the engagement from missile launch to the 
terminal flight phase where the missile is in the vicinity 
of the target is called the missile flyout. The Endgame 
phase includes the missile fuzing sequence for target 
detection and subsequent warhead detonation, followed by an 
evaluation of the effectiveness of the damage mechanisms 
associated with a warhead on the target under the mis- 
sile/aircraft encounter conditions at the time of warhead 
detonation . 

The overall measure of effectiveness of a missile system 
against an aircraft can be stated in terms of the single 
shot probability of kill, PKSS. The PKSS value encompasses 
all missile system functions from launch through the com- 
plete target engagement. significant factors considered in 
computing PKSS are launcher and missile reliability, gui- 
dance and control, fuzing, warhead detonation and evaluation 
of target damage and target kill. PKSS can be expressed in 
the general form as 

PKSS = Pi * Pg/1 * Pf/g * Pd/f * Ph/d * Pk/h 

where Pi is the probability of launch, Pg/1 is the probabil- 
ity (or occurance) ofsuccessful guidance given launch, Pf/g 
is the probability (or occurance) of fuzing given guidance, 
Ph/d is the probability (or number) of hits on the target 
given detonation, and ?k/h is the probability of target kill 
given the hits on the target. 

Missile Endgame programs measure the probability of 
kill, PK/E, of a terminal encounter. This is done by 
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assuming tha- the probability of missile launch and sucess- 
ful guidance (Pg/1 * Pi) is unity; therefore the missile is 
always in the vicinity of the target at the beginning of the 
encounter. The final result of the Endgame computations is 
an overall probability of kill given an encounter which can 
be expressed in the general form as 

PK/E = Pf/g * Pd/f * Ph/d * Pk/h 

The target detection, or fuzing, and the detonation are 
dependent upon the fuze operation and logic, the encounter 
geometry and the target signature. In many Endgame programs 
there are options to set Pf/g and Pd/f equal to unity in 
order to eliminate the fuze from consideration. 

The value of Ph/d, or the number and location of hits on 
the target, is computed by the program and is a function of 
the encounter geometry conditions, the target size and the 
missile warhead parameters. The Pk/h is determined by the 
design of the target model with it's associated component 
Pk/h functions or by the component vulnerable area (Av) 
tables for a given target, depending upon the type of End- 
game simulation used. These two values are utilized in com- 
puting Pk/d, the probability of kill given a detonation, 
which is of the general form 

PK/D = Ph/d * Pk/h 

The PK/E determined by the Endgame programs is not a 
single shot probability of kill PKSS. PKSS considers the 
entire engagement and hence involves many factors that are 
not accounted for in PK/E alone. The reliability of the 
missile and launcher, missile performance, the guidance and 
control, and many other factors are considered in computing 
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PKSS. It can therefore be seen that PK/E is only valid for 
the final terminal phase of an engagement and assumes all 
missile system functions prior to that point are perfect. 

The need for valid Endgame simulation is obvious. The 
cost of conducting extensive tests on every missile/aircraft 
combination is not only prohibitive, but very time consuming 
and may be impossible due to unavailability of the systems. 
The modeling of encounters does require extensive physical 
test data, but does not require entire new systems to be 
tested; only those components for which no data is currently 
available need be tested. The time and money required for 
one physical test can be used to produce literally thousands 
of simulation runs. Simulation can also be used in planning 
an actual test firing to optimize the probability of col- 
lecting the data desired. The use of Endgame simulation 
will allow designers to evaluate changes made to missile or 
aircraft systems much earlier in the design process, allow- 
ing the manufacturer to produce a product needing fewer 
modifications to meet both performance and survivability/ef- 
fectiveness specifications. A service interested in buying 
a new missile or aircraft will be able to conduct prelimi- 
nary "fly-offs'* between proposed designs prior to funding 
full scale prototypes in order to narrow the field of con- 
tenders to those meeting the desired specifications. The 
time required to go from preliminary design to an engineer- 
ing model can be reduced substantially if non-viable options 
can be reduced or eliminated through simulated encounters. 
The use of Endgame simulations to develop improved tactics 
and engagement methods for both missiles and aircraft could 
result in substantially increased effectiveness of forces 
presently available. New or improved enemy systems can be 
modeled quickly and, as additional data becomes available. 
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updated to provide new tactics to counter the threat. Sev- 
eral of the aforementioned uses of Endgame simulation are 
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II. EXISTING ENDGAME SIMULATION PROGRAMS 



A. GENERAL 

The description of each program presented in this Chap- 
ter will include discussion about the type of modeling used, 
the input format, program options, how PK is assessed, the 
program originator and the original intended use of the 
simulation. This chapter is intended to present the various 
programs in a brief summary form and no attempt will be made 
here to evaluate the relative merits of the individual End- 
game simulations. 

3. ATTACK 

The ATTACK program is the current version of the 
AIR-TO-AIR TERMINAL SIMULATION (NWC TN4565-1-70) [Ref. 1], a 
Naval Weapons Center, China Lake, revision of a methodology 
developed at the Naval Missile Center, Point Mugu. The 
documentation for ATTACK was published in June 1974. The 
objective of ATTACK, as stated in it’s User Manual, "is to 
predict the ability of a missile to detect and destroy an 
airborne target." To this end, the program provides a PK 
assessment for four damage mechanisms. These are (1) direct 
hits, (2) blast, (3) multiple fragment (structural) , and (4) 
single fragment (component) damage mechanisms. 

This program requires four target geometric representa- 
tions, one representation for each of the damage mechanisms. 
The fuzing portion of the program requires a fifth target 
representation which is dependent upon the fuze type 
selected. 

The program is intended to provide results for the fol- 
lowing purposes: 

(1) Weapon system evaluation 

(2) Warhead design 
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(3) Fuze optimization 

(4) Survivability studies 

(5) Trade-off studies 

The methodology for damage assessment is composed of the 
following classes: 

(1) Structural 

(a) direct hit model 

(b) blast model 

(c) multiple fragment model 

(2) Component 

(a) single fragment model 

The direct hit model consists of a target modeled using 
triangular plates (see Fig. 2-1) and a missile which is 
represented by a collection of points (see Fig. 2-2) . 

The missile trajectory is computed from a user specified 
encounter geometry. The program determines if one or more 
of the missile points will intersect the target and the time 
at which first contact will occur. If the first contact 
occurs before fuzing, a direct hit kill is assessed and 
kills by other possible damage mechanisms are not evaluated. 
If fuzing occurs prior to first contact then a '•preempted' 1 
direct hit is recorded and noted in the output, but other 
damage mechanism kills are evaluated. 

The blast model is composed of a group of cylinders and 
hemispher ically shaped end caps surrounding the target body 
and extremities (see Fig. 2-3) . The radius assigned to each 
of the blast cylinders is a function of both the strength of 
the included structure and the explosive charge of the war- 
head. The cylinder radii, which are determined external to 
the program in a seperate analysis, are scaled automatically 
to a user specified encounter altitude. If the warhead 
detonates within the volume of one of these blast cylinders, 
a blast kill occurs and no other damage mechanisms are 
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Fig. 2-1 ATTACK Direct Hit Model 



18 




Fig. 2-2 ATTACK Missile Representation 
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Fig. 2-3 ATTACK Blast Model 
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evaluated. If the warhead detonates outside the volume of 
the blast cylinder, no damage is credited to blast effects, 
and the program goes on to evaluate other damage modes. 

The warhead model in ATTACK uses the concept of fragment 
spray zones and fragment weight classes. The warhead may 
have up to ten polar zones with up to five fragment mass 
classes in each. The user may also define up to eight 
radial zones which allows simulation of nonsymmetric frag- 
ment sprays. For detonation of the warhead, the user can 
choose from eleven fuze logics (the number of fuze logics 
may vary from one installation to another) . A typical list- 
ing of fuze logics is given in Table II-1. 

The multiple fragment model for structural damage uses a 
segmented cylindrical target representation (see Fig. 2-4) . 
The program increments the centroid of the cylindrical seg- 
ments by the target velocity vector from the time of warhead 
detonation. The fragment ballistics are computed as a func- 
tion of: 



(1) 


Fragment 


ma ss 


(2) 


F rag ment 


shape 


(3) 


Fragment 


initial velocity 


(4) 


Fragment 


drag coefficient 


(5) 


Target range and aspect from warhead at detonation 


(6) 


Fragment 


and target flight paths 



The number of fragments with their associated energies 
which strike each cylindrical segment is determined by the 
location of the target segment within one or more of the 
warhead dynamic radial and polar zones. The energy density 
is computed and compared with a critical level of energy 
specified for that segment by the user. If the calculated 
value exceeds the specified energy density, a structural 
kill is assessed. 
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TABLE II-1 







LISTING OF FUZE LOGICS 


Logic 


1 


Semi-active doppler fuze 


Logic 


2 


Seai-active doppler fuze with 
signal stretcher 


Logic 


3 


Semi-active doppler fast track fuze 


Logic 


4 


Semi-aczive doppler fuze for intercept 
arm; fixed angle fuze 


Logic 


5 


Fixed angle active fuze 


Logic 


6 


IR fuze operation in pursuit mode 


Logic 


7 


Active fuze with fore and aft 
fixed angle fuze cones 


Logic 


3 


Passive fixed angle fuze 


Logic 


9 


$emi-active fuze with guard channel for 
intercept arm. Fixed angle for home 
on jam, fuze on jam 


Logic 


10 


Semi-active doppler with guard channel 
arm 


Logic 


1 1 


Instantaneous detection* 


* N ote : 


This 


option has been added to the NPS version 
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Fig. 2-4 ATTACK Structural Damage Model 
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The single fragment or component kill model consists of 
individual components, represented by spheres or points, 
located at appropriate positions with respect to the target 
coordinate system origin (see Fig. 2-5) . The computational 
process for pk used in this model is similar to that used in 
the multiple fragment case. The component (sphere) centroid 
location and radius are used to determine the fractional 
area (FRACT) of the component within a given polar and 
radial fragment spray zone. 

FRACT = area of the component hit by fragment spray (At) 
component presented area (Ap) 

The vulnerability of the components to the impact of frag- 
ments is measured by the component vulnerable area, Av. 
Vulnerable area tables are input for each component in the 
model as a function of fragment impact aspect angle, frag- 
ment mass and fragment impact velocity. The distance of the 
component centroid from the warhead origin at detonation 
(DIST)is used to determine exact fragment impact velocity 
and the striking azimuth and elevation angles for a specific 
weight class. These are used in conjunction with the vuln- 
erable area tables to compute the appropriate component 
vulnerable area (Av) . A fragment beam area (FA) within the 
polar and radial zone boundaries is computed at the dis- 
tance, DIST. The fragment spray density (RHO) is given by: 

a HO = Q/FA 

Where Q is the number of fragments in the weight class and 
polar/radial zone considered. The expected number of lethal 
hits (E) for the specified weight class is computed from: 

E = RHO * Av * FRACT 
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Component Locations 
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Fig. 2-5 ATTACK Single Fragment Vulnerable Components 




The expected number of lethal hits is accumulated for each 
polar zone, radial zone and fragment weight class. The com- 
ponent Pk is computed by the following equation: 

Pk (component) = 1.0 - exp (-S) 

The encounter geometry as shown in Figure 2-6 is speci- 
fied by the user. The missile may be oriented with respect 
to either the target or to a relative velocity vector. The 
user may either specify a missile miss distance or closest 
point of approach (CP A) or utilize the program to generate 
random miss distances from a Gaussian distribution. The 
user may also specify a standard deviation for the miss dis- 
tance distribution, and multiple trajectories may be simu- 
lated for a given scenario. 

C. SCAN 

SCAN is a digital computer program developed under the 
supervision of the Pacific Missile Test Center, Point Mugu 
[Ref. 2]. The documentation was completed in June 1976 
under the auspices of the Joint Technical Coordinating Group 
on Aircraft Survivability (JCTG/AS) . The objective of the 
SCAN Endgame simulation, as defined in the User Manual, is 
"to predict the probability that an aircraft will survive an 
attack by a missile armed with a warhead." A PK is computed 
for three cases: 

(1) Direct hit 

(2) Blast 

(3) Fragment damage 

This program can be used to provide data for: 

(1) Aircraft design 

(2) Aircraft survivabilty studies 

(3) Justification of new survivability features 

/ 
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tig, 2 6 ATTACK Encounter Geomet 



One major feature of this program is the complex geome- 
tric model of the target. In SCAN, the target model is com- 
posed of a series of components, where each component is 
represented by one or more boxes, polygons, or guadric sur- 
faces with bounding planes (e. g. finite length cylinders, 
ellipsoids, etc. linked together by logical .AND. /.OR. 
statements) . A sample model is shown in Figure 2-7. 

Sach component is assigned a Pk/n value based upon one 
of three types of vulnerability. The three types are: 

(1) Single fragment vulnerable 

(2) Energy density vulnerable 

(3) Area removal vulnerable 

For the first type, the measure of vulnerability is the 
probability of component kill given a hit by a fragment 
(Pk/h) . This is expressed as a constant term plus a linear 
function of fragment mass and of impact velocity in the 
form: 

Pk/h = PK ( 1 ) + P K (2) * M + PK ( 3) * V 

where P K ( 1 ) is a constant term 

P K ( 2 ) is the coefficient of mass 
PK(3) is the coefficient of velocity 
M is the fragment mass in grains 
V is the fragment velocity in feet per second 

For the second type of vulnerability, the energy density 
kill is expressed in terms of a required minimum area 
exposed to a threshold energy density level with a limiting 
fragment mass below which no computations are made. This 
type of kill probability is generally applicable to target 
structural members, whereas single fragment vulnerability is 
commonly used for components. 

For the third type of vulnerability, the measure of an 
area removed kill is defined by a minimum area removed, 
below which no damage occurs, and an area which, if removed. 
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Fig. 2-7 Geometric Representation of the Exterior of the F-16 Aircraft 



will cause complete failure. The kill probability is com- 
puted linearly between these two values. 

For each target component the user must specify a mater- 
ial type and thickness. The material type is chosen from a 
list of ten options provided by the program and specified in 
Table II-2. A component surface is designated as solid or 
hollow and as either an internal or external aircraft 
component. 

Each component’s vulnerability and susceptability type 
is chosen from a list of eleven options shown in Table II-3. 
It is also possible to define a component to be non-vulnera- 
ble to specific damage mechanisms. Specific components may 
be designated as infrared (I R ) sources allowing the simula- 
tion of IR fuzing mechanisms. 

The killing of an individual component may or may not 
cause a target kill. Consegua ntly , aircraft subsystems can 
be defined by linking components with logical .AMD. /.OR. 
statements, and aircraft systems can be composed of previ- 
ously defined subsystems. The components are identified by 
the order in which they were input for the geometric repre- 
sentation. This feature of the program can be used to 
define multiply vulnerable or redundant components. Various 
levels of target kill can then be specified in tarms of com- 
ponents, subsystems and systems. 

The SCAN blast model and warhead model are both similar 
to the ATTACK model. The SCAN fuzing model has only three 
options: 

(1) Instantaneous detection 

(2) Infrared (IR) fuzing 

(3) Single look-angle active fuze 

The program has three possible scenario choices availa- 
ble. In one, the user may define a trajectory by fixing the 
initial missile range from the target and the orientation of 
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TABLE II - 2 



LIST OF MAT ERIAL TYPES 



Material Input Cod 

Magnesium 10 

Aluminum 2024T 20 

Titanium Alloy 30 

Face Hardened Steel 40 

Mild Steel 50 

Hardened Steel 60 

Lexan 70 

Plexiglass 80 

Doron 90 

Bullet. Resistant Glass 100 
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TABLE II- 3 



VULNERABILITT/SUSCEPT ABILITY OPTIONS 



Option number 
1 
2 

3 

4 

5 



Option Description 

Energy density vulnerable 

Single fragment vulnerable 

Area removal vulnerable 

Nonvulnera ble to fragments, 
direct-hit vulneraole 

Nonvulnera ble to fragments 
and direct hit 
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IR source and nonvulnerable 
to fragments 



7 



Energy density vulnerable, 
invisible to EM fuze 



8 



Single frgament vulnerable, 
invisible to EM fuze 



9 



Area removal vulnerable, 
invisible to EM fuze 



10 



Nonvulnerable to fragments, 
invisible to EM fuze 



1 1 



Nonvulnerable, invisible to 
EM fuze 
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the warhead polar and radial zones into a number of elements 
containing fragments of the same class which are all tra- 
velling in approximately the same direction. A. representa- 
tive ray is generated to characterize the fragment of each 
element, and the motion of this characteristic fragment is 
simulated along a trajectory. This procedure can be very 
time consuming when the number of fragments is large or when 
the target is complex. In order to reduce the computational 
requirements, the user must provide limiting spatial parame- 
ters. These parameters are dependent upon the physical 
dimensions of the target. Limits are established at values 
which slightly exceed the target dimensions. Outside of the 
limiting values no fragment computations take place. 

SCAN has two graphic display programs available for use 
with the simulation. Several facilities, including the 
Naval Postgraduate School, have installed the capabilities. 
The 5PLGEN program is a preprocessor which will check all of 
the geometric target model information for proper type and 
limits, then display the target model on a Tektronics 4081 
graphics display terminal. The SPDRAW program uses fragment 
impact data generated by a SCAN simulation to display the 
hits on a displayed target model. 



D. SESTEM II 

SESTEM II is a digital computer program developed by the 
Aeronautical systems Division - Deputy for Development Plan- 
ning (ASD/XRHD) , Wright - Patterson Air Force Base [Ref. 3]. 
The documentation was published in May 1977. The objective 
of the SESTEM II Endgame simulation, as defined in the model 
description, is "to evaluate the terminal effectiveness of 
missiles with nonnuclear warheads against U.S. and foreign 
aerial targets." The terminal effectiveness is measured in 
terms of the terminal encounter single shot probability of 
kill ?K . The program is designed to evaluate missiles with 
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blast - fragment warheads. A PK is computed for three 
cases : 



(1) 


Direct hit 




(2) 


Fragment damage 




(3) 


Blast 




This program has been utiliz 


ed for: 


(D 


Preliminary warhead 


design and fuze optimization 


(2) 


Formulating tactics 


and countermeasures require- 



ments 



(3) Computerized air-to-air duel simulations 

(4) Reconstructing and analyzing selected Southeast 
Asian combat incidents 

(5) Evaluating existing and conceptual aircraft in sup- 
port of an Offensive Air Support Mission Analysis 

The program requires three general types of input data: 

(1) Encounter data 

(a) terminal geometry (see Fig. 2-8) 

(b) missile aimpoint 

(c) target and missile encounter altitude 

(2) Missile warhead and Fuze data 

(a) circular error probable (CEP) , miss 
distance, or iso-Pk conditions 

(b) fuzing equations 

(c) fuze delay time 

(d) fragment sprayband and fragment density 

(e) fragment average mass and initial velocity 

(f) fragment cross-sectional area and 
coefficient (af drag 

(3) Target data 

(a) component size and location 

(b) individual component fragment vulnerable areas 

(c) external blast kill contours 
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Fig. 2-8 SESTEM II Terminal Geometry 



The target is simulated in the computer as a collection 
of shapes, to be discussed later, representing fragment 
vulnerable, masking, or fuzing components. Each fragment 
vulnerable component has assigned it's own appropiate table 
of fragment vulnerable areas as a function of aspect angle 
and impact velocity just as in the ATTACK program. The pro- 
gram constructs an external blast kill contour for the tar- 
get, missile and altitude being evaluated using the input 
data. The program represents each component by means of a 
grid of variable side length inscribed on the component sur- 
face. A "target point" is generated in the center of each 
grid square, represented by direction cosines and X, Y, Z 
coordinates. The target point is then used to represent 
that grid square in fragment interaction computations. The 
missile warhead parameters and fuzing equations are simu- 
lated in the program using the static input data. Dynamic 
resolution of the static warhead data is done by the pro- 
gram. Various types of fuzes, e.g. radar, contact, proxim- 
ity, may be simulated. Both the target and the missile are 
assumed to be flying constant speed, straight-line trajecto- 
ries during the terminal phase. The missile and target 
approach each other along their relative trajectories until 
the fuzing equations are satisfied. After the appropiate 
delay time, warhead detonation occurs and the dynamic inter- 
action with the vulnerable components is computed. The 
probabilities of killing the target by blast, direct hit, 
and each component by fragments, are computed and combined 
to predict the probability of target kill. Fuzing effects 
may be examined by selecting as many as eleven fuzing points 
for each trajectory. 

Trajectories may be generated by various methods with 
S2STEM II. Parallel trajectories may be generated randomly 
by assuming a bivariant normal or other type of expected 
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distributions. Sequential groups of up to thirty trajecto- 
ries may be generated and averaged. Discrete trajectories 
may also be input and the PK computed for each case. Aver- 
age PK values as a function of miss distance can be calcu- 
lated by specifying various miss distances. The program 
locates trajectories on a circle of radius equal to the spe- 
cified miss distance, calculates the individual trajectory 
PK values and the average PK for this miss distance. 

The program can also be used to generate "iso-PK" cont- 
ours. In this mode the program computes the blast kill 
boundary (PK =1), the no fuzing boundary (PK = 0) and PK 
for varying distances between these two values. 

The target is modeled using ellipses and parallelo- 
pipeds. A fuselage or fuselage-shaped component with a 
cross section parallel to the target Y-Z plane is repre- 
sented by the upper and lower halves of two ellipses with 
common minor axes. The center of the ellipses may be dis- 
placed an arbitrary distance along a line parallel to the 
target Z-axis. Each component may be represented by using 
up to thirty cross sections (see Fig. 2-9). The equations 
defining the surface of the component and the direction 
cosines of any target point on the surface are generated by 
t.he program as required. A wing or wing-like component with 
a cross-section parallel to the X- Z plane of the target is 
represented by ellipses as shown in Figure 2-10. The verti- 
cal stabilizer or similarly configured components with 
cross-sections parallel to the X-Y plane are represented by 
ellipses also. The program will compute the direction 
cosines for a point on the surface of these structures as 
well as the equations for the surfaces generated. The 
program has the capability to represent components such as 
fuel tanks or electronics as rectangular parallelopipeds . 

The parallelopipeds are located by specifying the corner 
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Fig. 2-9 SESTEM II Fuselage Model 
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Fig. 2-10 SESTEM II Wing Model 
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with the lowest target coordinate values and X, Y, Z dis- 
placement values. The program has the capability to gener- 
ate "mirror image" components by reversing the Y-values. 

This allows easy representation of axisymmetric components 
using only one side as input. 

The expanding blast wave from the warhead detonation can 
cause target destruction from the overpressure and dynamic 
effects on the structure. A blast contour encloses that 
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Fig. 2-11 SESTEM II Blast Contour 




2-12 SESTEM II Warhead Model 
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SHAZAM utilizes component Pk/h values to compute target 
vulnerability to fragments and uses a ray trace method to 
determine fragment paths from detonation. 

The SHAZAM program uses a target model of a fixed for- 
mat. The exterior surface of the target is a finite element 
model composed of polygons. This allows SHAZAM to utilize 
the target models generated by the FASTGEN or SHOTGEN compu- 
ter programs. The interior components are modeled using 
lines, spheres, cylinders, and polygons. 

Blast can be handled statically or dynamically by 
SHAZAM. Each external polygon can have a seperate blast 
kill radius specified to accurately model blast vulnerabil- 
ity. The program is capable of adjusting the shape of the 
blast contours with time for dynamic blast modeling. 

The program can combine the blast, fragment, and direct 
hit Pk values to obtain an overall target pk. 

SHAZAM has the capability to generate graphics and uti- 
lize interactive graphics. This is an added option that is 
not currently part of the "core" program. 

F. HEFMOD 

The REFMOD digital computer program is a reference model 
used for computing the effectiveness of externally detonat- 
ing weapons against moving targets. The model was developed 
under the auspices of the Joint Technical Group for Muni- 
tions Effectiveness (JTCG/ME) , Anti-Air Missile Evaluation 
group. The first version, REFMOD-1, had documentation com- 
pleted in December 1979. The latest version, REFMOD-2, has 
preliminary documentation dated March 1981 [Ref. 4], 

REFMOD-2 will be the version described in this section. 

REFMOD has been assembled by incorporating methodologies 
from other existing Endgame simulations, such as AMEGS, 
ATTACK, SCAN, SHAZAM, and WHDEVAL. In order to combine fea- 
tures from these programs, it was necessary to modify them 
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to provide consistant nomenclature and coordinate systems 
for the resultant program. 

REFMOD also includes some significant additional fea- 
tures which enable it to work with a wide variety of vulner- 
ability models and allows it to evaluate warhead/target 
combinations that were previously too cumbersome to assess. 
The fuzing routine utilized will allow greater flexibility 
in fuze modeling. 

The program is currently being rewritten in Standard 
FORTRAN using a structured programming format. This will 
result in a program usable at any facility with Standard 
FORTRAN capability and enough computing power. 

The REFMOD program has a high degree of flexability. 
Three different modes of encounter input are available. The 
fuzing can either be performed by the program or determined 
externally, and any convential warhead type can be modeled. 
Three shape options are available for contact hit target 
modeling. There are three options for blast kill computa- 
tions, and several types of fragment damage can be modeled. 

The types of studies that REFMOD can be used for are: 

(1) To evaluate the operational effectiveness of exist- 
ing and proposed missile systems. 

(2) To assist in fuze optimization and selection of 
warhead design. 

(3) To assess aircraft survivability and countermeasure 
effectiveness. 

The type of target model required as input data depends 
upon the damage mechanism specified in the vulnerability 
model. The direct hit model utilizes the physical shape of 
the target to determine if the missile contacts the target. 
The target model can be described by any combination of the 
three shapes listed: 

(1) Truncated elliptical cone 
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(2) Polygon 

(3) Ellipsoid 

Figure 2-13 gives examples of direct hit modeling. A full 
range of target Pk values from zero to one can be specified 
for a contact or debris kill. 

The evaluation of blast effects can use three different 
model types for evaluation. The three models are: 

(1) The target blast model is represented by hemispher- 
ically capped cylinders. A kill is recorded according to 
whether the warhead detonation occurs inside at least one of 
the blast cylinders (see Fig. 2-14) . 

(2) The target is modeled using blast ellipsoids and a 
kill is recorded according to whether the warhead detonation 
occurs inside at least one of the blast ellipsoids. 

(3) The target is modeled by blast ellipsoids and blast 
damage centers. Data on blast pressure and impulse are 
required to determine blast kills. This model provides for 
ncn-spherical blast representation. 

The area removal model for fragment kills consists of a 
representation of individual structural members considered 
vulnerable to fragmentation effects. The probability of 
removal of a specified amount of material in a "contigious" 
path across a structural member is determined based upon 
fragment striking conditions, fragment spray angles, and 
other pertinent information. This value is used to deter- 
mine the probability of structural member failure and is 
combined with the probabilities of failure of all other 
structural members to give a total kill probability. Figure 
2-15 is an example of this model. 

There are four other fragment vulnerable component types 
that can be used. Each component can use only one type for 
vulnerability modeling. The four models are: 

(1) Vulnerable area 



46 



A. THUNCATEO ELLIPTICAL CONES 





Fig. 2-13 REFMOD Direct Hit Modeling 
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Fig. 2-14 REFMOD Blast Cylinder Model 
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Fig. 2-15 REFMOD Structural Failure Model 
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(2) Function of mass/velocity /density 

(3) Expected value 

(4) Focused Fragment Controlled Motion (FFCM) 

For the vulnerable area type, components can be 

described as being spherical, linear, cylindrical, or planar 
in shape (see Fig. 2-16) . Each of these shapes is used with 
the traditional vulnerable area methodology, using the vuln- 
erable area tables generated by the C0VA3T computer program. 

For the mass/velocity/density type the components are 
described as cylinders (see Fig 2-17) . The component kill 
condition is given by: 

F = A * (M ** B) * (V ** C) per unit area 

Where A, B, and C are constants to be defined by the user, K 
is fragment mass, and 7 is fragment impact velocity. This 
vulnerability measure can vary as a function of the fragment 
impact angle, and is described by upper and lower threshold 
boundaries as illustrated in Figure 2-18. The conditional 
kill probability is zero when F is below the lower threshold 
value. When F is above the upper threshold level, the kill 
probability is set equal to one. The component PK is inter- 
polated linearly between the two threshold values. Energy 
density (A = 0.5, B = 1, C = 2) and momentum density (A = 1, 
3 = 1, C = 1) are specific examples of this type of vulner- 
ability model. 

The expected value model represents the components as 
cylinders or line segments that outline the components. The 
vulnerability varies as a function of impact angle. The 
vulnerability is described in terms of Pk/h for N fragments, 
where H varies from one to a maximum specified value. The 
number of hits is calculated using a Poission distribution. 
This model type can also be used to specify a lethal radius 
vulnerability criterion. 
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Fig. 2-16(a) REFMOD Spherical Vulnerable Component 
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Fig. 2-16(b) REFMOD Linear Vulnerable Component 
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Fig. 2-16(c) REFMOD Cylindrical Vulnerable Component 
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Fig. 2-16 (d) REFMOD Planar Vulnerable Component 
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Fig. 2-17 REFMOD Mass/Velocity /Density Vulnerable Model 
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Fig. 2-18 REFMOD Mass/Velocity/Density Function 



The Focused Fragment Controlled Motion (FFCM) model is 
the same as the expected value model except that the 
fragments are assumed to be uniformly spaced instead of uni- 
formly distributed. The uniform spacing of the fragments is 
due to the use of large preformed fragments with predictable 
ejection angles. These preformed fragments maintain the 
same angular spacing with respect to the origin of the war- 
head detonation throughout their flight. 

The REFMOD program is capable of handling the following 
warhead types: 

(1) Continuous rods - the ends of the rods are welded 
together causing the rods to open in a hoop. 

(2) Divergent fragmentation - most wide beam warheads 
are of this type. 

(3) Convergent fragmentation - multiple point initia- 
tion causes the fragment patterns to cross over one another. 

(4) Focused Fragment Controlled Motion - large pre- 
formed fragments with highly predictable ejection angles, 
resulting in a uniformly spaced pattern instead of a uni- 
formly distributed one. 

(5) Aimable - special cases of the above types in which 
the fragment density is non-uniform about the roll axis. 
Figure 2-19 illustrates the warhead types described above. 

Several fuze routines have been revised and documented 
for use with REFMOD. These routines simulate the fuzes cur- 
rently employed on existing missile systems. REFMOD is 
compatible with fuzing data obtained from one of these 
routines, from routines designed for conceptual fuzing, from 
fuzing data obtained from flight tests, or from any other 
source. The fuzing routine has a special target model made 
up of line segments representing surfaces that can be sensed 
by the fuze (see Fig. 2-20). 
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Fig. 2-19 REFMOD Warhead Types 
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Fig. 2-20 REFMOD Fuzing Model 
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III. COMPARISON OF FEATURES AND TECHNIQUES 



This chapter will compare the various techniques used in 
computer Endgame simulations. 

A. GEOMETRIC MODELING; FINITE ELEMENT VERSOS COMBINATORIAL 

The two types of target modeling commonly use in Endgame 
programs are finite element (polygons) or combinatorial 
geometry (COMGEOM) . 

Finite element modeling uses multi-sided polygons to 
construct the target and/or its components (see Fig. 3-1). 
The modeling accuracy can vary greatly, depending upon the 
types and number of polygons used or allowed. An advantage 
to this type of modeling is that there are computer programs 
available that are designed to generate this type of target 
model (e.g. FASTGEN or SHOTGEN) . All of the programs dis- 
cussed in this study use a finite element model for comput- 
ing the results of at least one of their damage modes. 

CCMGEOM modeling involves the use of not only polygons, 
but also boxes, cylinders, hypreboloids, ellipsoids, and 
various conic sections with bounding planes. These complex 
shapes allow very exact modeling of a target and its compo- 
nents (see Fig. 3-2) . However, the use of these complex 
geometric shapes adds a degree of difficulty to accurate 
target modeling. It appears that there are currently no 
computer programs designed to generate COMGEOM target 
models. One advantage to most programs that utilize COMGEOM 
models is that they can also accept finite element modeling 
as input. The SCAN program utilizes COMGEOM target models. 
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Fig. 3-1 Finite Element Model 
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C0MGE0M Model 



and both SESTEM II and SHAZAM utilize some of tha features 
of COMGEOM target modeling. 

B. VULNERABLE AREA VERSOS COMPONENT PROBABILITY OF KILL 

The two techniques used in assessing fragment damage are 
vulnerable area and component probability of kill given a 
hit, Pk/h. 

1 . Vulnerable Area 

The vulnerable area method utilizes vulnerable area 
tables for each fragment vulnerable component. Each compo- 
nent will have as many vulnerable area tables as there are 
fragment mass classes to be considered. The standard type 
of vulnerable area table has 26 entries per fragment mass 
class per component for eight velocities, based upon the 26 
aspect angles shown in Figure 3-3 (for example, ATTACK 
requires 26 input data cards per component for one fragment 
mass class) . 

Components that are not fragment vulnerable, but act 
as masking or shielding for other components, must also have 
a full set of null vulnerable area tables in order to ascer- 
tain the amount of fragment slowdown or deflection. 

The use of vulnerable area modeling of components 
simplifies the amount of computing required during program 
execution. The program computes the aspect angle of the 
fragments from the user input data, and then uses the input 
value of fragment mass to find the correct vulnerable area 
values from the tables input by the user. The programs usu- 
ally use linear interpolation for aspect angles that fall 
between those tabulated. ATTACK, REFMOD and SESTEM II all 
use vulnerable area modeling. 

2 . Component Probability of Kill Given a Hit 

The method of using component probability of kill 
given a hit, Pk/h, involves assigning each component modeled 
a set of Pk/h values. The values assigned to the component 
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Fig. 3-3 Aspect Angles for Vulnerable Area 
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are usually in a functional relationship and are used to 
determine component kill levels. For example, SCAM uses a 
function of the form 

Pk/h = PK ( 1 ) + P K ( 2) * H ♦ PK ( 3) * V 

where P K ( 1 ) is a constant term 

?K(2) is the coefficient of mass 
PK(3) is the coefficient of velocity 
M is the fragment mass in grains 
V is the fragment velocity in feet per second 

The modeling required for component Pk/h is often 
complex. The components can be modeled with a high degree 
of accuracy to obtain valid results. The shapes used in 
modeling the target vary depending upon the program being 
utilized, but the SCAN program discussed earlier is a good 
example of how many possible geometric shapes can be used. 
The Pk/h assigned to a component has no aspect dependency 
that has to be tabulated as in vulnerable area modeling. 

The component Pk/h type of Endgame program does all the com- 
puting of aspect angle corrections internally. 

A positive feature of the complex model is that it 
can serve a multitude of functions. The same geometric 
model can be used to assess other damage types such as 
direct hit, structural, and blast damage. 

This method of modeling gives more flexabiiity to 
the user. The Pk/h functional relationship for component 
kill assessment can be easily changed and the modification 
of individual components does not require an outside program 
or creation of a new input table of data. 

Programs using component ?k/h modeling are capable 
of evaluating the effect of spallation. Since each com- 
ponent can have a seperate material type and thickness spe- 
cified, it is possible to generate fragments of various 
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materials , sizes, and velocities from a single warhead 
fragment. 

The amount of computing time required for this type 
of modeling is often greater than that required for a vuln- 
erable area model; however, by limiting the spatial volume 
being considered by the program to the immediate vicinity of 
the target (thus not computing the paths of non-str iking 
fragments) , the computing time can be reduced substantially. 

The difficulty level of constructing targets manu- 
ally can vary greatly depending upon the target complexity 
and the geometric shapes used in modeling. 

C. WARHEADS AND FUZING 

All of the warhead models assume definable polar and 
radial zones of fragment ejection. The number of zones used 
and the number of fragment mass classes that may be con- 
tained within any one zone differs in the various Endgame 
programs. Only REPMOD permits the modeling of continuous 
rod and Focused Fragment controlled Motion (FFCM) warheads. 
The SCAN warhead model is the only one that allows fragments 
to be composed of various material types and permits a 
choice of fragment shapes. 

The programs considered here utilize two general methods 
to determine the fuze type. SCAN, ATTACK and SESTEM II all 
have option lists from which the user choses a specific fuze 
type. REF MOD and SHAZAM allow the user to model any fuzing 
desired in a subroutine and utilize that subroutine in the 
program to fuze the warhead. 

Vulnerable area type Endgame programs often require a 
seperate model for the fuze initiation (see Fig 3-4). Com- 
ponent ?k/h type Endgames usually have the fuzing model 
incorporated into the geometric model. 
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Fig. 3-4 Vulnerable Area Fuze Model 
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D. INPOT FORMAT 

The ease of inputting information to the program can 
greatly affect the overall "usability" of an Endgame simula- 
tion. This section will discuss several methods of input- 
ting the necessary information. 

All of the programs use some type of formatted input. 

The input data is supposed to be in a predefined order, with 
specific values placed in designated columns on an I3M card 
(or a facsimile file on disk or tape storage) . 

One method of input is to use one data file for all of 
the information required. This requires the manipulation of 
a large block of data in order to make minor changes. It 
also requires the user to be familiar with the entire input 
file to make changes and assemble data for program execu- 
tion. ATTACK, REFMOD and SESTEM II all utilize this type of 
input. 

Another method is to subdivide the input data into 
smaller, functionally oriented input files. The SCAN pro- 
gram, as it is implemented at the Naval Postgraduate School, 
is an example of this technique. The function oriented 
files deal with a specific part of the input required such 
as warhead data, encounter geometry, and the target geome- 
tric model. This allows the user to modify a specific sec- 
tion of input without having to manipulate all of the input 
data. It also allows the compilation of a library of func- 
tional data files that may be combined to provide the 
desired scenario for a simulated encounter. This method is 
more user oriented, and while the programming required for 
this type of input may be more complex (not always the 
case) , the ease of use justifies the time spent in addi- 
tional programming. 

One major feature that should be in any Endgame program 
is a preprocessing graphics capability for plotting the 
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targets . Since it is very easy to make mistakes when pre- 
paring the target modeling data, especially in the COMGEOM 
approach, a plot of the modelled target is essential for 
verification. Only SCAN currently has this capability. The 
SPLGEN program takes the same input for the geometric model 
as SCAN does. It uses the input data to check for correct 
modeling and then produces a plot or display of the model. 

3. AVAILABILITY OF TARGET MODELS AND OTHER INPOT DATA 

The two programs that have the most documented input 
data and target models available are SCAN and REFMOD . SCAN 
is a component Pk/h type program utilizing a COMGEOM target 
model and REFMOD is a vulnerable area type program using a 
finite element direct-hit model. Since both models have 
been extensively used, large quantities of both types of 
data are available. New and updated targets and data are 
continuously being produced for use with both SCAN and 
REFMOD. The ATTACK program can utilize the REFMOD data with 
only minor modifications. Information about the production 
of data and models for the other programs was not readily 
available. 

F. PROGRAM OQTPUT 

The output format must be user oriented if the program 
is to be useful. The computer must be utilized to do as 
much of the output organization and interpretation as 
possible. This does not mean it should supress any of the 
output data, but rather the computer should be used to pro- 
cess and collate the raw data in such a manner as to allow 
direct use of the output without having to do additional 
interpretation. All of the programs considered in this 
study incorporate this ideal to varying degrees. The abil- 
ity of these programs to allow changes in output format var- 
ies from one simulation to the next. 
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The ability to combine damage types and create systems 
or subsystems from the target component output is an essen- 
tial reguirement. All of the programs have the capability 
to do this, but the SCAN program is the easiest to use. 

SCAN allows the linking of individual components into sys- 
tems through the use of logical .AND/. OR statements. 



G. GRAPHICS CAPABILITIES 

The need to generate graphic output from Endgame simula- 
tions has been sorely underemphasized. A preprocessor 
graphics program can be invaluable to the target modeler in 
tracing errors in the input data file, and the added impact 
of a visual representation of the fragment hit locations 
using a post- processor is immense. The ability to view the 
damaged area of the target and to follow the fragments as 
they pass through the target can be a tremendous aid when 
interpreting the output data. 

Only the SCAN and SHAZAM programs contain a graphics 
option. The version of SCAN used at the Naval Postgraduate 
School is linked to two graphics programs, a preprocessor 
and a postprocessor. The use of graphics at N?s allowed the 
correction of several complex target model inaccuracies that 
never would have been discovered without the added capabili- 
ties. Figure 3-5 is an example of computer generated 
graphics. 

H. COMPUTATIONAL TIME AND COST REQUIREMENTS 

It should be obvious that the faster the program exe- 
cutes, the better it seems from the standpoint of the user. 
The multitude of factors that control program execution time 
are, in most cases, dependent upon the computer facilities 
available and not upon the programs themselves. Since at a 
large modern computer facility, the execution times of even 
the larger simulations are under five seconds (single 
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Fig. 3-5 Computer Generated Graphics 



encounter case) , the question of execution time is really 
more a question of the cost of computer time. The rela- 
tively high cost of computer time at many facilities could 
make a vulnerable area type of program more appealing due to 
a slightly faster execution time (this would be dependent 
upon the specific case ) . The use of any graphics programs 
can also be very time consuming and therefore costly. Since 
the cost-per-computati on of computer time is dropping due to 
advances in technology, the cost differential between diffe- 
rent types of programs should decrease. 

I. DOCUMENTATION 

With the exception of SCAN, the documentation for the 
current Endgame programs is wholly inadequate. Most of the 
programs had example problems to use for familiarization, 
but in several cases the example values given did not match 
the parameters of the User Manual. In all cases the docu- 
mentation for the math models of the programs was very good. 
The key weaknesses were failure to explain how the actual 
program functioned and to thoroughly explain the input data 
parameters . 

The SCAN documentation suffered from very few of these 
problems. It is user oriented and easily understood. In 
addition to an overall example problem, SCAN also has many 
small examples that illustrate individual parameters. There 
is a complete explanation of how the program functions and 
how each subroutine works. 

The use of explanatory comments in the body of the com- 
puter program itself is another problem area. All of the 
programs had some explanation in the program text, but in 
all cases it was not sufficient to allow a new user to 
quickly locate and understand individual program functions. 
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IV. RECOMMENDATIONS FOR FUTURE ENDGAME PROGRAMS 



This chapter will discuss those features that seem best 
suited to a standardized, user oriented Endgame simulation. 

A. BASIC REQUIREMENTS 

The programming should be done in a standard language 
that is compatible with all types of computers. The most 
logical choice would be the most current version of Standard 
FORTRAN. This will allow the program to be machine indepen- 
dent and avoids the present problem of program translation. 

The program decided upon should be for use by all ser- 
vices and not oriented toward any special user group. This 
will eliminate the current tendency towards "private" End- 
game programs that cannot readily be compared to one 
another. It would also simplify the problems cf corpora- 
tions that are working or bidding on service contracts. All 
companies would have the same Endgame simulation so that 
comparisons between bidders on the basis of simulation 
results would be valid, and a company or research facility 
could use the same simulation regardless of the particular 
service they were doing work for. 

A single tri-service program would also generate a vast 
data base that could be referenced by any qualified user. 
This would avoid the current problem of constructing multi- 
ple, independent data bases that are not transferable due to 
program incompatabil ity . This would require the creation of 
an organization to manage and maintain the program and it’s 
data base. 

The use of a single program also means that improvements 
would be disseminated quickly and that the program should 
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converge rapidly to a more valid simulation due to many 
facilities working with and trying to optimize the same 
program. 

B. TARGET MODELING AND VULNERABILITY TYPE 

The use of combinatorial geometry (COMGEOM) modeling is 
the most efficient choice. By using the COMGEOM model of 
the target for fragment, energy density, area removal, and 
direct hit vulnerability, only one additional model, the 
blast model, would be required. If the blast radius was 
modeled for each component as it is in the SHAZAM program, 
only one model would be required for all damage modes. 

The selection of COMGEOM would still allow the use of 
finite element modeling and the associated computer gener- 
ated targets, but it would also be able to use the more 
accurate complex geometric shapes. Since it is feasible 
that a program could be written to generate COMGEOM target 
models, it could eventually be possible to have both the 
accuracy of COMGEOM and the finite element ease of modeling. 
Until such a program is developed, a master library of 
models and components would eliminate most of the need for 
independent model development. 

The use of component Pk/h for fragment vulnerability is 
well suited for use with COMGEOM target modeling since in 
many cases an entire component can be modeled with one 
shape. The use of component Pk/h also allows the user to 
modify the Pk/h of individual components, possibly to 
reflect new physical test data. This will also permit the 
accurate modeling of spallation effects by allowing the use 
of various material types in constructing the target model, 
and permit modification of the function used to determine 
component Pk/h. Additionally, the COMGEOM component Pk/h 
model also permits the use of many material types in 
simulating components, an important factor in view of the 
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rapid changes in aircraft structural materials presently 
occuring. 

C. WARHEAD AND FOZE MODELING 

The warhead model should be able to handle the numerous 
types of warheads in use or proposed for use. The most 
efficient way to handle the warhead model is not by using 
built-in warhead options, but instead by using specialized 
warhead subroutines that are input with the other data and 
called by the main program as is done in SHAZAM. This would 
allow very accurate warhead models to be created and permit 
new warhead designs to be added at a later date. It would 
streamline the main program by not having to offer an exten- 
sive selection of warhead options. 

The fuzing models should also be subroutines input with 
the data for the same reasons that apply to the warhead 
model . 

D. INPUT FORMAT 

The input should be broken up into small functional sub- 
groups, e.g. fuzing, warhead, target model etc. This type 
of input allows easy manipulation of data and simplfies 
tracing errors in input data. 

E. OUTPUT FORMAT 

The output format utilized by SCAN is a good choice to 
use. It mirrors the input values and encounter conditions, 
lists user defined systems and subsystems, tabulates PK 
information for each component and the user defined systems, 
and gives a summary table for each type of kill mechanism 
for the target as a whole. 
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P. FLYOUT SIMULATION INTERFACE 

The ultimate missile versus aircraft program would be a 
missile flyout simulation that interfaced with an Endgame 
simulation. This interfacing of simulations would allow the 
computation of PKSS values for the missile systems simu- 
lated. SHAZAM has been used in this manner so that this is 
not an untried concept. Ideally, the flyout and Endgame 
programs should be designed as one large program that encom- 
passes two "stand alone" programs. A truely sucessful inte- 
gration of this type would require the generation of flyout 
programs that generate output specifically structured to 
provide the values required by the Endgame program. 

G. GRAPHICS CAPABILITY 

It is essential that new Endgame simulations provide for 
graphic output, preferably ineractive graphic routines. The 
two graphic routines designed for use with SCAN, SPLGEN and 
SPDRAW, are good examples of the graphics programs needed by 
Endgame simulations. One of these routines, SPLGEN, is a 
preprocessor that checks the target model for geometric 
accuracy than plots the model on a plotter or a cathode ray 
tube (CRT) display. The preprocessor is independent of the 
actual Endgame program, and is extremely useful to someone 
modeling targets for Endgame simulation use. The other 
routine, SPDRAW, is a postprocessor that utilizes the output 
from the Endgame program. The Endgame program generates a 
data file containing the location of all the fragment hits 
on the target model. This routine plots the aircraft then 
superimposes the fragment hits onto the aircraft plot. This 
program also displays output on either a plotter or a CRT 
display . 

The visual image of the target damage is extremely use- 
ful in interpreting the damage caused by a missile warhead. 
The use of graphics will allow extensive investigation of 
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Fig. 4 - 1 ( a ) Exterior View Graphics Example 
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Fig 0 4-1 (b) Interior Components Graphics Example 
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Fig. 4-l(c) Exterior View Graphics Showing Hits 



It is important that any changes made to the program be 
thoroughly documented and overall program documentation 
updated on a regular basis. If the program is not kept veil 
documented as changes occur, the whole purpose of documenta- 
tion is defeated and an unfamiliar user will not be able to 
utilize the program without having to do extensive trial and 
error testing to find the undocumented changes. 



81 



V 



DIRECTION OF CURRENT RESEARCH 



This chapter will discuss those programs that are still 
in development or are undergoing major modification. This 
information was provided by the program developers. 

A. SCAN 

The SCAN program is currently in the process of being 
modified at the Pacific Missile Test center, Point Mugu. 

The new version is referred to as "Son of SCAN" and has been 
substantially streamlined. Many of the modifications being 
incorporated are to enhance the graphics capability of the 
original program. It is not known at this time whether Son 
of SCAN will become generally available. 

3. SHAZAM 

The SHAZAM program, by intent, will always be in a state 
of modification. The documentation package for SHAZAM 
should become available within the next year. The documen- 
tation will explain the "bookkeeping" core of the program 
and discuss the requirements for subroutines to be added by 
the user. 

C. REFMOD 

The latest version of REFMOD, REFMOD II, is not fully 
operational at this time. Modifications are being carried 
out at the Naval Weapons Center, China Lake. The latest 
version has a more flexible fuze package, more warhead 
options, better documentation, and a much simpler input for- 
mat. Preliminary documentation is currently available, and 
the new program should be available by August 1981. 
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APPENDIX A 



A USER'S MANUAL FOR SCAN AT NPS 

Tnis User's Manual provides instruction for the execu- 
tion of SCAN at NPS using the IBM 3278 display terminal. 

A. STORAGE REQUIREMENTS 

The compiled version of the SCAN program requires 
approximately one cylinder of storage on the I3M 3033 system 
(2240 records) . The source code version (FORTRAN) requires 
approximately two cylinders of storage (4322 records). 
Execution of the program typically requires less than one 
cylinder of storage for the output. The exact amount of 
storage required during execution is dependent upon the num- 
ber of cases, the number of fragments in the warhead, and 
the complexity of the target 

3. USER INSTRUCTIONS FOR EXECUTION 

This section will deal with the execution of SCAN using 
existing case, warhead, and target geometry files stored on 
a read-only disk file. The creation of new input files will 
be discussed in later sections. 

1. Turn on the terminal. 

2. When the large "NPS" logo appears on the screen, 
press the RESET key followed by the ENTER key. 

3. When "CP READ" appears in the lower right of the 
screen, type in L nnnnP, where nnnn is your user identifica- 
tion number. Press ENTER. 

4. You will now be asked for your password. Type in 
your password (it will not appear on the screen), then press 
ENTER. 
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5. Type "CP LINK TO xxxxP 191 AS 192 HR", where xxxx 
is the user number of the project file. Press ENTER. 

6. You will now be asked for the project password. 

Type it in (it will not appear on the screen) . Press ENTER. 

7. Type "ACCESS 192 B" and press ENTER. 

8. Type "PROFILE EXEC" and press ENTER. This step 
assumes that you do not currently have a PROFILE EXEC file 
on your private disk. If you have your own PROFILE EXEC, 
you must ensure that it contains the following command: 

GLOBAL TXTLIB FORTMOD2 M0D2EEH 

9. The SCAN program will require three data files in 
order to execute properly. These files are: 

(a) GEOM DATA - a file containing the geometric 
target model, limiting parameters and kill expressions. 

(b) WARHEAD DATA - a file containing the missile, 
warhead and fuzing parameters. 

(c) CASE DATA - a file containing the encounter 
geometry information. 

The WARHEAD DATA and CASE DATA files are normally prese- 
lected and ready for use. Changes to these files will be 
discussed in later sections. The GEOM DATA file must be 
chosen from a list of available targets (see Table AA-1) . 

To transfer the target file chosen to the GEOM DATA file for 
input into the program, type in "COPY Fn Fm 31 GEOM DATA 
Al"and press ENTER. Fn refers to the file name of the tar- 
get desired and Fm refers to the file mode of the target 
file as given in Table AA-1. 

10. You are now ready to run SCAN. At NPS, the SCAN 
program is stored under the title "SCANMAIN". This version 
has been specially modified for the NPS computer system. 

You should insure that there is enough room on your disk for 
the output files (approximately one cylinder) . 
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TABLE AA-1 



TARGETS AVAILABLE FOR USE WITH SCAN 



Target 


Filename 


file 


A-7 


A7DATA 




Harpoon 


A3M8 6 


DATA 


Tomahawk 


A3 M 1 09 


DATA 


3ackf ire 


BACK F 


DATA 


Shoebox 


BOX 


DATA 


Drone 


BQM3 4 


DATA 


Drone 


B2M1 07 


DATA 


Exoce t 


EXOCET 


DATA 


Foxb at 


FOXBAT 


DATA 


Kingf ish 


KINGF 


DATA 


Kitchen 


KITCH 


DATA 
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You 



12. Type in "RUN SC ANMAIN " then press ENTER. 



should 


get 


the : 


following listing on your screen: 




FILEDEF 


01 


DISK 


SCAN 1 FORTRAN 


(RECFM FA BLOCK 


131 


PERM 


FILEDEF 


05 


DISK 


CASE DATA 








FILEDEF 


06 


DISK 


SCAN 2 FORTRAN 


(RECFM FA BLOCK 


131 


PERM 


FILEDEF 


11 


DISK 


GEOM DATA 








FILEDEF 


12 


DISK 


WARHEAD DATA 








FILEDEF 


19 


DISK 


IMPACT DATA (RECFM V3S LRECL 


127 


BLOCK 131 


FILEDEF 


30 


DISK 


FILE DATA 









LOAD SCANMAIN 
START 

EXECUTION BEGINS... 

The program will take a few seconds to execute. The actual 
run time will depend upon computer workload, target and war- 
head complexity, and the number of cases per run. 

12. Once the program has completed running, type "L" 
and press ENTER to see if the output files were generated on 
your disk. The output files that should be created are: 

(a) SCAN1 FORTRAN A1 - this file is an echo print 
of target parameters. 

(b) SCAN2 FORTRAN A1 - this file contains a sum- 
mary of the encounter conditions, warhead/missile parame- 
ters, system definition statements and output summary. 

(c) FILE DATA Al - this file contains the system 
definition statements in assembly language. It is created 
and used by the SCAN program, and is not of any value to the 
user. 
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(d) IMPACT DATA A1 - this file contains the assem- 
bly language ceding for the fragment impact points. It is 
used in the SCAN graphics program. For information concern- 
ing the use of graphics with the SCAN program, see the the- 
sis of LCDS T. Hayes [Ref. 5]. 

13. To obtain printouts of the SCAN results, type 
"PRINT SCAN1 FORTRAN (CC" then press ENT SR then type "PRINT 
SCAN2 FORTRAN (CC" and press ENTER. The printout of these 
two files will be output by the line printer in Ingersoll 
140 and filed alphabetically in the output bins according to 
your last name. Figures A-1 and A-2 are examples of SCANl 
FORTRAN A 1 and SCAN2 FORTRAN A1 printouts. 

14. To view the output of the program at the terminal 
display screen, use the XEDIT mode of the NPS VM/CMS time 
sharing system to review the files created (a guide to the 
full capabilities of XEDIT is available in the consultants 
office at the computer center) . It will be necessary to 
move the display right or left to view all of the informa- 
tion in the files due to the line length of the output (the 
method for doing this is explained in the XEDIT User's 
Guide) . 

C. INSTRUCTIONS FOR MODIFYING OP. CREATING A WARHEAD 

This section deals with the creation of new WARHEAD DATA 
files and the modification of existing WARHEAD files for use 
with the SCAN Endgame program. For modeling purposes, the 
CG of the target is assumed to be at the origin of the tar- 
get coordinate system as shown in Figure A-3. In reality, 
this may not be the actual target CG. Check the target 
model to find the origin of the target coordinate system by 
utilizing the target plots available in the thesis of LCDR 
T. Hayes [Ref. 5] before you start in order to avoid input- 
ting incorrect data values. 
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Fig. A- 1 SCAN1 FORTRAN Example output 
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Fig. A - 2 ( a ) SCAN 2 FORTRAN Example Outp 



vivu Nil tvirnhWH-i n**iH'iVP 






— a 
' r -v 




t * 




90 



Fig. A - 2 ( b ) SC AN 2 FORTRAN Example Output 
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Fig. A- 2(c) SCAN2 FORTRAN Example Output 
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Fig. A-3 Assumed Target CG 
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6. In columns 1 - 10: enter the number of static polar 
zones in your warhead (up to 36). This is an integer value. 

7. In columns 11 - 20: enter the number of fragment 
mass classes for each static polar zone (up to 3) . This is 
an integer value. Press ENTER. You will note that the 
first line has moved up and the cursor is now ready for the 
next line. 
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8. In columns 1 - 10: enter the lower angle of the 
first static polar zone (see Fig. A-4) for definitions of 
lower and upper angles). A real value of F10.3 format, 
angle measured in degrees from the forward end of missile 
roll axis (0 - 180) . 

9. In columns 11 - 20: enter the upper angle of the 
first static polar zone. A real value, F10.3 format, mea- 
sured in degrees from the forward end of missile roll axis 
(0 - 180) . 

10. In columns 21 - 30: enter the speed of the first 
mass class of fragments at the lower boundary of the first 
static polar zone. A real value, F10.3 format, in units of 
feet/sec. 

11. In columns 31 - 40: enter the speed of the first 
mass class of fragments at the upper boundary of the first 
static polar zone. A real value, F1Q.3 format, in units of 
feet/sec. 

12. In columns 41 - 50: enter the mass of the first 
mass class of fragments ejected in the first static polar 
zone. A real value, F10.3 format, in units of grains. 

13. In columns 51 - 60: enter the total number of frag- 
ments of the first mass class contained in the first ststic 
polar zone. A real value, F10.3 format. 

14. In columns 61 - 70: enter the initial position of 
the center of fragments of the first mass class in the first 
static polar zone with respect to the center of the warhead, 
as measured along the missile roll axis. A real value, 

F10.3 format, the units are feet. 

15. In columns 71 -74: enter the material code indicat- 
ing the type of material for the first mass class of frag- 
ment in the first static polar zone. Table AA-2 lists the 
material types available and the appropriate code for each 
type. An integer value, 14 format. 
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Fig. A-4 Warhead Static Polar Zones 
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TABLE AA-2 



LIST OF MATERIAL TYRES 



Material Input Cod 

Magnesium 10 

Aluminum 2024T 20 

Titanium Alloy 30 

Face Hardened Steel 40 

Mild Steel 50 

Hardened Steel 60 

Lexan 70 

Plexiglass 80 

Doron 90 

3ullet Resistant Glass 100 
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16. In columns 77 -80: enter the shape of the fragments 
in the first mass class in the first static polar rone. 

Table A- 3 lists the shapes available for use. An alphanu- 
meric character string, A4 format. Press ENTER. 

17. Repeat steps 8-16 for each fragment mass class in 

the first static polar zone then repeat steps 8-16 for 
each additional static polar zone (up to 36), and repeat 
steps 8-16 again for each fragment mass class in each 
static polar zone. The number of times you will have to 
execute steps 8 - 16 is equal to the number of fragment mass 
classes multiplied by the number of static polar zones. 
EXAMPLE: for three fragment mass classes and two static 

polar zones, steps 8-16 will be repeated six times. 

18. In columns 1 - 10: enter the type of fuze to be 
simulated. The options are: 

(a) 0 for instantaneous detection and detonation. 

(b) 1 for fuze on IR source only (requires IR 
source on target) . 

(c) 2 for fuze on any reflected target source. 

An integer value, 110 format. 

19. In columns 11 - 20: enter the position of the prox- 
imity fuze target detection device (TDD) with respect to the 
warhead center. This is a real value, F10.3 format, mea- 
sured in feet along the missile roll axis. 

20. In columns 21 - 30: enter the delay time between 
target detection and warhead detonation. A real value, 

F10.3 format, measured in seconds. 

21. In columns 31 - 40: enter the mean value of the 

proximity fuze cone half-angle (look-angle) as measured from 

the forward missile roll axis. A real value, F 10.3 format, 

measured in degrees (0 - 180) . 

22. In columns 41 - 50: enter the standard deviation of 

the fuze cone half-angle, assuming a normal distribution of 
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TABLE AA-3 



FRAGMENT SHAPES AVAILABLE FOR USE WITH SCAN WARHEAD MODEL 



Fr 



agment Shape 
Cube 

Spheroid 
Rectang ular 
Irregular 



Program Code 
CUBE 
SPHE 
RECT 
IRRE 
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angles. A real value, F10.3 format, measured in degrees 
(0 - 180) . 

23. In columns 51 - 60: enter the proximity fuze cutoff 
range. A real value, F10.3 format, measured in feet along 
the normal to the roll axis from the TDD position. If the 
fuze type selected was 0 then enter 0.0 for this value. 

24. In columns 61 - 66: enter the radius of the missile 
cylindrical body. A real value, F6.2 format, measured in 
feet. 

25. In columns 67 - 72: enter the position of the mis- 
sile contact fuze (usually at the tip of the nose) . A real 
value, F6.2 format, measured in feet from warhead center 
along the missile roll axis. 

26. In columns 73 - 78: enter the distance the missile 
extends aft of the warhead center. A real value, FS.2 for- 
mat, measured in feet along the missile roll axis. Press 
2NT2R . 

27. In columns 1 - 10: enter the fuselage blast radius. 
This is the maximum distance from the target centerline at 
which detonation of the warhead will cause catastrophic 
structural failure, assuming the target is at sea level. A 
real value, F10.2 format, measured in feet from the target 
centerline . 

28. In columns 11 - 20: enter the distance from target 
CG to the front of the fuselage blast cylinder. A real 
value, F10.2 format, measured in feet. 

29. In columns 21 - 30: enter the distance from target 
CG to the back of the fuselage blast cylinder. A real 
value, F10.2 format, measured in feet. 

30. In columns 31 - 40: enter the wing blast radius. 
This is the maximum distance from the wing at which 
detonation of the warhead will cause catastrophic structural 
failure, assuming the target is at sea level. A real value, 
F10.2 format, measured in feet from the wing centerline. 



99 



31. In columns 41 - 70: snter the X, Y, Z components of 
the end point of the wing blast centerline closest to the 
target fuselage. The X, Y, Z values are real, 3F10.2 for- 
mat, measured in feet from the target CG. Press ENTER. 

32. In columns 1 - 30: enter the X, Y, Z components of 
the end point of the wing blast centerline furthest from the 
target fuselage. The X, Y, Z values are real, 3F10.2 for- 
mat, measured in feet from the target CG . Press ENTER. An 
illustration of the blast model is shown in Figure A-S. 

33. You have now completed your warhead design. Press 
ENTER. The cursor should now be back in the lower left cor- 
ner of the screen and an END OF FILE statement should show 
up as the last entry in the file you just finished. Type 
'•FILE" and press ENTER. Your screen should now revert to 
the format it had before you entered the XEDIT mode. Type 
"L" and ENTER. You should see a listing of all the files on 
your disk space. The warhead DATA file you created should 
now appear in that listing. 

34. To use your warhead design in the SCANMAIN program, 
the warhead file must be named "W ARHEAD" . If the file you 
created is not named WARHEAD, you must rename it. Type in 
"COPY ?n DATA A1 WARHEAD = =" then press ENTER. The Fn 
refers to your original filename. This will not rename your 
original file, instead it creates a copy of it with the cor- 
rect name. Thus, on your file listing you will have your 
original warhead file and a duplicate of it under the file- 
name of WARHEAD. 

D. MODIFYING A WARHEAD FILE 

This section will deal with modifying an existing 
WARHEAD DATA file. Extensive reference will be made to Sec- 
tion C of this Manual. 

1. Turn on the terminal and log onto the system. 
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Fig . 



A-5 SCAN Blast Model 
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2. Type "XEDIT Fn Fm" and press ENTER. Fn refers to 

the filename of the file you wish to modify, and Fm refers 
to the filemode, usually a DATA file. An example is: XEDIT 

WARHEAD 1 DATA. 

3. You should now have the data in the file visible on 
the display screen. If the data covers more than one 
screen, use the following commands to survey the file: 

(a) Press "ALT'* and "PF8" keys at the same time to 
advance one screen forward in the file. 

(b) Press "ALT" and "PF8 " keys at the same time to 
revert back one screen in the file. 

4. Use Section C of this Appendix to locate the posi- 
tion of the values you wish to change. Display the section 
of the data file to be changed on the screen by using 
instruction 3 above. 

5. By using the four cursor positioning keys just to 
the right of the main keyboard, position the cursor under 
the values you wish to change. Simply type the new values 
in over the previous ones. Make sure that the new values 
are justified correctly and in the right format as discussed 
in Section C. 

6. After completing your changes on the page on the 
screen, press ENTER. The cursor will return to the lower 
left corner of the screen and the changes will have been 
made. You must enter all changes on each full screen of 
data before "paging" through the file. 

7. Once all changes have been entered and the cursor 
is back in the lower left corner of the screen, type "FILE" 
and press ENTER. This will permanently file your change and 
return you to the normal operating mode of the system. The 
corrected file is now ready for use. 
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3. CREATING A NEW CASE DATA FILE 

This section will deal with the creation of new CASE 
DATA files for use with the SCAN Endgame simulation. The 
information in this section is also applicable to the modi- 
fication of CASE DATA files currently available as discussed 
in the next section. 

1. Turn on the terminal and log onto the system. 

2. Type in "XEDIT xxxxxxxx DATA” and press ENTER, 

x xxxxxxx is the filename (up to eight spaces, the first 
space must be a letter, numbers may be used in the others) . 
An example is: XEDIT TESTCASE DATA. 

3. You are now in the XEDIT mode of operation. A 
guide to the full capabilities of XEDIT is available from 
the consultants office at the computer center. 

4. Type "I” and press ENTER. This will put you in the 
input mode required to create a new file. 

5. You will see an index line across the center of the 
screen. The index numbers correlate to the columns on an 
IBM data card. The input line is directly below the index 
line which will help you place the required data in the cor- 
rect columns for input. 

***** ** *************** ******** ************** ************** 

IMPORTANT NOTE!!! 

All integer values must be RIGHT JUSTIFIED in the field spe- 
cified. Real values may be anywhere in the specified field. 
Letter character strings must be RIGHT JUSTIFIED in the 
field specified. 

* **#5* ** ## jjC ****** ## A*#*#*#**#*****#*#**#********#**** 

6. In columns 1 - 10: enter the type of missile tra- 
jectory desired from Table AA-4. An integer value, 110 for- 
mat. Press ENTER. 
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TABLE AA-4 



MISSILE TRAJECTORY OPTIONS 



Cods Meaning 



0 


Terminate execution of program 


1 


Indicates a fixed trajectory or 
detonation point specified by an 
initial position measured from 
the target C3 


2 


Indicates a trajectory with a fixed 
missile guidance error (CPA) 


3 


Indicates trajectory in which the 

missile guidance error is 
computed from a normally 
distributed sample with a given CEP 
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7. In columns 1 - 10: enter the number of missile tra- 
jectories to be considered. An integer value, 110 format. 

8. In columns 11 - 20: enter the target speed at time 
of intercept. A real value, F10.3 format, measured in 
feet/second. The value must be greater than zero. 

The steps 9, 10, and 11 refer to Figure A-6. 

9. In columns 21 - 30: enter the roll angle (PSI) of 
the target at intercept. A real value, F10.3 format, mea- 
sured in degrees. The roll angle is measured with respect 
to the horizontal "flat earth" reference plane, right wing 
"down" is positive (0 - 360) . 

10. In columns 31 - 40: enter the pitch angle (THETA) 
of the target at intercept measured with respect to the flat 
earth plane. A real value, FI 0.3 format, measured in 
degrees. A positive angle indicates a climb, a negative 
angle, a dive (-90 to +90) . 

11. In columns 41 - 50: enter the yaw or heading angle 
(PHI) of the target at intercept. A real value, F10.3 for- 
mat, measured in degrees from the Y-axis of the reference 
plane. positively increasing yaw is measured in a counter- 
clockwise rotation (0 - 360) as viewed looking down the 
Z-axis of the reference coordinate system toward the origin. 

12. In columns 51 - 60: enter the missile speed at 
intercept. A real value, F10.3 format, measured in feet/ 
second (>0) . 

13. In columns 61 - 70: enter a mean value for missile 
angle of attack. A real value, F10.3 format, measured in 
degrees. Figure A-7 depicts the angle desired. 

14. In columns 71 - 80: enter the standard deviation of 
the missile angle of attack, assuming a normal distribution 
of angles. A real value, F10. 3 format, measured in degrees. 
Press ENTER. 
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Fig. A-6 Target Roll, Pitch and Yaw Angles 
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Missile Angle of Attack 



15. In columns 1 - 10: enter the mean pitch angle of 
the missile. A real value, FI 0.3 format, measured in 
degrees (-90 to +90) with respect to a "flat earth". A 
positive angle indicates a climb and a negative angle, a 
dive (see Fig. A-8) . 

16. In columns 11 - 20: enter the standard deviation of 
the missile elevation angle, assuming a normal distribution 
of elevation angles. A real value, F10.3 format, measured 
in degrees. 

17. In columns 21 - 30: enter the mean value of the 
azimuth angle of the missile in the terminal intercept. A 
real value, F10.3 format, measured in degrees (0 - 360) and 
referenced to the target coordinate system if target roll, 
pitch, and yaw are set at zero. If the target has roll, 
pitch, and yaw values the azimuth angle of the missile is 
referenced to the flat earth reference coordinate system 
(see Fig. A-9) . If the roll, pitch, and yaw of the target 
are at zero, a value of zero for missile azimuth angle 
implys a tail-chase encounter, a value of 180 implys a 
head-on encounter. 

18. In columns 31 - 40: enter the standard deviation of 
the missile azimuth angle assuming that the distribution of 
angles is normal. A real value, F10.3 format, measured in 
degrees . 

19. In columns 41 - 50: enter the altitude above sea 
level at which the encounter takes place. A real value, 
F10.3 format, measured in feet. 

20. In columns 51 - 80: enter the X, X, Z components of 
the missile aimpoint with respect to the target CG. A set 
of real values, 3F10.3 format, measured in feet. This is 
the point on the target the missile is aiming for. Miss 
distances in the output are calculated with respect to the 
aimpoint selected. Press ENTER. 



108 




Fig. A-8 Missile Pitch and Azimuth Angles 
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Fig. A-9 SCAN Coordinate System Relationships 
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For user specified trajectories or point detonation of 
the missile warhead, do steps 21-25. 

For simulations to determine the average survival prob- 
ability for a fixed guidance miss distance, do steps 26-28. 

For situations in which the user desires the simulation 
to generate the initial engagement geometry from a distribu- 
tion of encounter conditions and in which the miss distance 
for individual trajectories is drawn from a bi-variant nor- 
mal distribution of specified CEP, do steps 29-31. 

21. In columns 1 - 30: enter the X, Y, Z components of 
the initial position of the missile (for instantaneous deto- 
nation this is the point where the warhead will explode) 
measured in the target coordinate system (with respect to 
the target CG) . A set of real values, 3F10.3 format, mea- 
sured in feet. 

22. In columns 31 - 40: enter the target angle of 
attack with respect to the target velocity vector. A real 
value, F10.3 format, measured in degrees, counterclockwise 
positive. 

23. In columns 41 - 50: enter the target sideslip angle 
with respect to the target velocity vector. A real value, 
F10.3 format, measured in degrees, the counterclockwise 
direction is positive. 

24. In columns 51 - 60: enter the missile angle of 

attack with respect to the missile velocity vector. A real 

value, F10.3 format, measured in degrees, counterclock- 
wise positive. 

25. In columns 61 - 70: enter the missile sideslip 
angle with respect to the missile velocity vector. A real 
value, F10.3 format, measured in degrees, positive counter- 
clockwise. Press ENTER -and go to step 32. 

26. In columns 1 - 10: enter the closest point of 

approach of the missile trajectory to the specified aimpoint 
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on the target. A real value, F10.3 format, measured in feet 
from the missile aimpoint radially outward. 

27. In columns 11 - 20: enter the target angle of 
attack with respect to the target velocity vector. A real 
value, F10.3 format, measured in degrees, counterclockwise 
positive. 

23. In columns 21 - 30: enter the target sideslip angle 
with respect to the target velocity vector. A real value, 
F10.3 format, measured in degrees, positive counterclock- 
wise. Press ENTER and go to step 32. 

29. In columns 1 - 10: enter the circular error proba- 
ble (CEP) of the missile. A real value, F10.3 format, mea- 
sured in feet. The CEP will be centered about the missile 
aimpoint by the program. 

30. In columns 11 - 20: enter the target angle of 
attack with respect to the target velocity vector. A real 
value, F10.3 format, measured in degrees, counterclockwise 
positive. 

31. In columns 21 - 30: enter the target sideslip angle 
with respect to the target velocity vector. A real value, 
F10.3 format, measured in degrees, positive counterclock- 
wise. Press ENTER. 

32. You have now completed your CASE DATA file. press 
ENTER. The cursor should now be in the lower left corner of 
the screen and an END OF FILE statement should show up as 
the last entry in the file you just finished. Type "FILE” 
and press ENTER. Your screen should now rever to the format 
it had prior to entering the XEDIT mode. Type "L" and press 
ENTER. You should now see a listing of all the files on 
your alloted disk space. The file you just created should 
now be included in that listing. 
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F. MODIFY ING A CASE FILE 

This section will deal with modifying an existing CASE 
DATA file. Extensive reference will be made to Section E of 
this Manual. 

1. Turn on the terminal and log onto the system. 

2. Type "XEDIT Fn Fm" and press ENTER . Fn refers to 

the filename of the file you wish to modify. Fm refers to 
the filemode of the file you wish to modify (usually DATA). 
An example is: XEDIT CASE1 DATA. 

3. You should now have the data in the file visible on 
the display screen. If the data covers more than one full 
screen, use the following commands to survey the file: 

(a) Press "ALT" and "PF8" keys at the same time to 
advance one screen forward in the data file. 

(b) Press "ALT" and "PF8" keys at the same time to 
revert back one screen in the data file. 

4. Use Section S of this Manual to locate the posi- 
tions of the values you wish to change. Display the section 
of the data file to be changed using instruction 3 above. 

5. By using the four cursor positioning keys just to 
the right of the main keyboard, position the cursor under 
the values you wish to modify. Simply type the new values 
in over the previous ones. Make sure that the new values 
are justified correctly and in the right format as shown in 
Section C. 

6. After completing your changes on the page on the 
screen, press ENTER. The cursor will return to the lower 
left corner of the screen and the changes will have been 
made. You must enter all changes on each full screen of 
data before "paging" through the file. 

7. Once all changes have been entered and the cursor 
is back in the lower left corner of the screen, type "FILE" 
and press ENTER. This will return you to the normal 
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operating mode of the system and store the corrected file 
back on your disk space. The corrected file is now ready 
for use. 

G. OTHER MODIFICATIONS AND EXAMPLES 

Modifications to the GEOM DATA files, which contain the 
target geometric model, limiting parameters and kill expres- 
sions, are also possible. Changes to these parameters are 
often very complex and are not advised for someone unfamilar 
with the mechanics of the SCAN program. Information con- 
cerning these changes can be found in the User Manual for 
the SCAN program (master copy, not modified for NPS) . 

Figure A- 1 0 shows typical WARHEAD DATA and CASE DATA 
files. 

Figure A— 1 1 is a quick reference guide to the WARHEAD 
file input values. 

Figure A-12 is a quick reference guide to the CASE file 
input values. 
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10 WARHEAD DATA and CASE DATA Examples 
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Fig. A-ll(a) WARHEAD DATA Quick Reference Guide 
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Fig. A-ll(b) WARHEAD DATA Quick Reference Guide 
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Fig. A-12(b) CASE DATA Quick Reference Guide 



APPENDIX B 



PROGRAM CHANSES X APE TO SCAN 

This Appendix contains the changes made to the SCAN End- 
game program at the Naval Postgraduate School in order to 
function on the IBM 3033 computer system and to output data 
for the SP DRAW computer graphics program. several changes 
were also made in order to simplify the task of interpreting 
the output. The changes made will be listed and referenced 
by the routine or subroutine in which they occur. 

A. MAIN PROGRAM 

1. Added the following: 

REWIND 19 

This ensures that the data file used for graphics is stored 
correctly. 

B. BLOCK DATA 

1. Changed the DATA NAMARR array from: 

DATA NAMARR/4HMAGN, 4HSIUM, 4 HA LOM, 4 H2 024 , 4HTIIN, 4H ALOY, 

1 4HSTEE,4HFACE, 4HSTES, 4HM ILD, 4HSTEE,4 HHARD,4 HLEX A , 1H , 

2 4HPLX-,4HGLAS, 4HDORO,1H , 4H3ULL, 4HRES / 
to: 

DATA NAMARR/4HMA GN , 4 HSIO M, 4HALUM, 4H2024 , 4 HTITN, 4 HALOY, 

1 4HFACE,4 HSTEL, 4H MI LD, 4 HS TEL, 4HH ARD , 4 H ST EL , 4 HLEX N , 1 H , 

2 4HPLXI, 4HGLAS, 4HDORO, 1H ,4HBTJLL, 4HRESG/ 

2. Changed the DATA 70LTYP array from: 

DATA 7ULT YP/4HEN ER , 4 HDEN S , 4HS F ,4H7UL , 4 H AREA , 4 HR MV L, 

1 4HNON , 4HVCJL ,4HNV ,4HDHIT,4HNV ,4HIR S,4HED , 

1 4 HRRTP, 4 HS F , 4HRRTP, 4 HA R , 4HRRTP ,4HN7 ,4HRRTP, 

1 4HDHIT,4 HRRTP/ 
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t o : 

DATA VULTYP/4HENER, 4HDENS, 4 HSI NG , 4 H F RAG , 4 HAREA , 4 HRMV L, 

1 4 H , 4 HNON ,4HNONV / 4HDHII,4HNONV / 4HIR S,4HED R, 

1 4 HTR NS, 4 HSF R, 4HTRNS,4HA R R , 4HTRNS , 4 K N V R ,4 HTRNS , 

1 4HNVDH, 4HTRNS/ 

3. Changed the DATA SOT YP array from: 

DATA SKNTYP/4HEXTE,4HSKIN, 4 HEXTE, 4 HS LI D , 4 HINTE, 4 HSKIN, 

1 4 HINT E, 4 HSLID/ 

to : 

DATA S KNTYP/4HEXTE ,4HSKIN,4HEXTE,4HS0LI,4HINTE,4HSKIN, 

1 4HIHTE,4HS0LI/ 

4. Changed DATA TYP(2) from: 

DATA TYP ( 2) /4HHLPD/ 

to: 

DATA TYP (2) /4HELIP/ 

5. Changed DATA TYP(4) from: 

DATA TYP (4) /4HHC0N/ 

to: 

DATA T YP (4) /4HELC0/ 

C. SUBROUTINE READIN 

1. Added the following: 

COMMON/FCTR/FACTOR 

This is a value used to correctly dimensionalize the output 
data for graphics. 

D. SUBROUTINE UPDATE 

1. Added the following: 

COMMON/FCTR/FACTOR 

This is a value recieved from Subroutine READIN to correctly 
dimensionalize the output data for graphics. 

2. Added the following: 

REAL RANC, TPC 
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These are output values. 

3. Changed the following: 

DIMENSION R A (3) , RANDHT (3) , RANM (3) / VBAR(3) , VX(3) 
to : 

DIMENSION R A (3) , RANDHT (3) , RANM (3) ,VBAR(3) ,VX(3) , 

RANC (3) , TPC ( 3) 

4. Added the following: 

INTEGER I HIT 

An output value. 

5. Added the following: 

REAL FRAGQ,TFRST,THIT 

These are output values. 

6. Added the following: 

TFRST = 0. 

THIT = 0. 

IHIT = 0 

FRAGQ = 0. 

This initializes the new variables. 

7. Added the following: 

R ANC ( 1) = RAN ( 1 ) / FACTOR 

R ANC ( 2) = RAN ( 2) / FACTOR 

RANC (3) = RAN (3) / FACTOR 

TPC ( 1) = TP ( 1) / FACTOR 

TPC (2) = TP (2) / FACTOR 

TPC (3) = TP (3) / FACTOR 

WRITE (19) FRAGM, TFRST, RANC,TH IT, IPC/IHIT 

This outputs the fragment mass and it's coordinates as it 
enters a component. 

8. Added the following: 

RANC ( 1 ) = RAN ( 1 ) / FACTOR 

RANC (2 ) = RAN (2) / FACTOR 



122 



RANC (3) = RAN (3) / FACTOR 

TPC (1) = TP (1) / FACTOR 
TPC (2) = TP (2) / FACTOR 
TPC (3) = T P ( 3) / FACTOR 

WRITS (19) FRAGH ,TFRST, R AN C, T H IT , TPC , IH IT 

This outputs the fragment mass and it's coordinates as it 
exits a component. 

9. Added the following: 

FRAGQ = -999. 

RANC ( 1) = RAN (1) / FACTOR 

RANC (2) = RAN (2) / FACTOR 

RANC (3) = RAN (3) / FACTOR 

TPC (1) = TP (1) / FACTOR 

TPC (2) = TP (2) / FACTOR 

TPC (3) = TP (3) / FACTOR 

WRITS (19) FRAG2/TFRST,RANC,THIT,T?C,IHIT 

This outputs the end-of-file data for the graphics program. 

S. FUNCTION ASIN (X) 

1. Changed FUNCTION ASIN (X) from: 

FUNCTION ASIN (X) 

ASIN = ARSIN(X) 

RETURN 
END 
to : 

FUNCTION AFNSN(X) 

AFNSN = ARSIN(X) 

RETURN 

END 

This was required for compatibility with the I3M 3033 compu- 
ter system. All calls for ASIN(X) in the entire program 
have been changed to AFNSN (X) . 

F. File definition statements for SCAN 
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1. The file definitions in effect when using SCAN are: 

FILEDEF 01 DISK SCAN 1 FORTRAN 

FILEDEF 05 DISK CASE DATA 

FILEDEF 06 DISK SCAN2 FORTRAN 

FILEDEF 11 DISK GEOM DATA 

FILEDEF 12 DISK WARHEAD DATA 

FILEDEF 19 DISK IMPACT DATA 

FILEDEF 30 DISK FILS DATA 

These file definitions are contained in an executive program 

labeled SCANMAIN EXEC. 
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